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Abstract 

Let P be a set of n points in ]R * I 2 3 , and let DT(P) denote its Euclidean Delaunay triangulation. 
We introduce the notion of an edge of DT(P) being stable. Defined in terms of a parameter 
a > 0, a Delaunay edge pq is called a-stable, if the (equal) angles at which p and q see the 
corresponding Voronoi edge epq are at least a. A subgraph G of DT(P) is called {col, oc)-stable 
Delaunay graph (SDG in short), for some constant c > 1, if every edge in G is a-stable and every 
ca-stable of DT(P) is in G. 

We show that if an edge is stable in the Euclidean Delaunay triangulation of P, then it is 
also a stable edge, though for a different value of a, in the Delaunay triangulation of P under 
any convex distance function that is sufficiently close to the Euclidean norm, and vice-versa. In 
particular, a 6a-stable edge in DT(P) is a-stable in the Delaunay triangulation under the distance 
function induced by a regular /c-gon for k > In/a, and vice-versa. Exploiting this relationship 
and the analysis in f3|, we present a linear-size kinetic data structure (KDS) for maintaining 
an (8a, a)-SDG as the points of P move. If the points move along algebraic trajectories of 
bounded degree, the KDS processes nearly quadratic events during the motion, each of which 
can processed in O(logn) time. Finally, we show that a number of useful properties of DT(P) 
are retained by SDG of P. 
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1 Introduction 

Let P be a set of n points in IR 2 . For a point p E P, the (Euclidean) Voronoi cell of p is defined as 

Vor(p) = {x E R 2 | ||xp|| < Hop'll Vp' E P}. 

The Voronoi cells of points in P are nonempty, have pairwise-disjoint interiors, and partition the 
plane. The planar subdivision induced by these Voronoi cells is referred to as the (Euclidean) 
Voronoi diagram of P and we denote it as VD(P). The Delaunay graph of P is the dual graph of VD(P), 
i.e., pi] is an edge of the Delaunay graph if and only if Vor(p) and Vor(i]) share an edge. This is 
equivalent to the existence of a circle passing through p and q that does not contain any other 
point of P in its interior—any circle centered at a point of 9Vor(p) n 9Vor (q) and passing through 
p and q is such a circle. If no four points of P are cocircular, then the planar subdivision induced 
by the Delaunay graph is a triangulation of the convex hull of P—the well-known (Euclidean) 
Delaunay triangulation of P, denoted as DT(P). See Figure [l] (a). DT(P) consists of all triangles 
whose circumcircles do not contain points of P in their interior. Delaunay triangulations and 
Voronoi diagrams are fundamental to much of computational geometry and its applications. See 
HI for a very recent textbook on these structures. 




Figure 1 . (a) The Euclidean Voronoi diagram (dotted) and Delaunay triangulation (solid); (b) the Q -Voronoi diagram and 
Q-Delaunay triangulation for an axis-parallel square Q, i.e., the Voronoi and Delaunay diagrams under the Loo-metric. 


In many applications of Delaunay/Voronoi methods (e.g., mesh generation and kinetic collision 
detection), the input points are moving continuously, so they need to be efficiently updated 
as motion occurs. Even though the motion of the points is continuous, the combinatorial and 
topological structure of VD(P) and DT(P) change only at discrete times when certain "critical 
events" occur. A challenging open question in combinatorial geometry is to bound the number of 
critical events if each point of P moves along an algebraic trajectory of constant degree. 

Guibas et al. Ifl4j showed a roughly cubic upper bound of 0(n 2 A s (n)) on the number of critical 
events. Elere A s (n) is the maximum length of an (n, s)-Davenport-Schinzel sequence 12H . and s is a 
constant depending on the degree of the motion of the points. See also Fu and Lee Ifl3l . The best 
known lower bound is quadratic |2T| . Recent works of Rubin [19, 20] establish an almost quadratic 
bound of 0 (h 2+£ ), for any e > 0, for the restricted cases where any four points of P can be cocircular 
at most two or three times. In particular, the latter study Il20l covers the case of points moving along 
lines at common unit speed, which has been highlighted as a major open problem in discrete and 
computational geometry; see Ifl2l . Nevertheless, no sub-cubic upper bound is known for more 
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general motions, including the case where the points of P are moving along lines at non-uniform, 
albeit fixed, speeds. It is worth mentioning that the analysis in Ifl9l . and even more so in Il20l . is 
fairly involved, which results in a huge implicit constant of proportionality. 

Given this gap in the bound on the number of critical events, it is natural to ask whether one 
can define a large subgraph of the Delaunay graph of P so that (i) it provably experiences at most a 
nearly quadratic number of critical events, (ii) it is reasonably easy to define and maintain, and 
(iii) it retains useful properties for further applications. This paper defines such a subgraph of the 
Delaunay graph, shows that it can be maintained efficiently, and proves that it preserves a number 
of useful properties of DT(P). 

Related work. It is well known that DT(P) can be maintained efficiently using the so-called 
kinetic data structure framework proposed by Basch et al. IflOl . A triangulation 7 of the convex 
hull conv(P) of P is the Delaunay triangulation of P if and only if for every edge pq adjacent to 
two triangles A pqr + and Apqr in 7, the circumcircle of A pqr + (resp., A pqr~) does not contain 
r (resp., r + ). Equivalently, 

Zpr + q + Zpr~q < n. (1) 

Equality occurs when p, q, r + , r~ are cocircular, which generally signifies that a combinatorial 
change in DT(P) (a so-called edge flip) is about to take place. We also extend {1} to apply to edges 
pq of the hull, each having only one adjacent triangle, A pqr + . In this case we take r to lie at 
infinity, and put Zpr~q = 0. An equality in |Tjl occurs when p, q, r 1 become collinear (along the 
hull boundary), and again this signifies a combinatorial change in DT(P). 

This makes the maintenance of DT(P) under point motion quite simple: an update is necessary 
only when the empty circumcircle condition (HI fails for one of the edges, i.e., for an edge pq, 
adjacent to triangles A pqr + and A pqr~, p, q, r , and r~ become cocircularj’J Whenever such an 
event happens, the edge pq is flipped with r 1 r ~ to restore Delaunayhood. Keeping track of these 
cocircularity events is straightforward, and each such event is detected and processed in 0(log n) 
time Ml4l . However, as mentioned above, the best known upper bound on the number of events 
processed by this KDS (which is the number of topological changes in DT(P) during the motion), 
assuming that the points of P are moving along algebraic trajectories of bounded degree, is near 
cubic Il4l (except for the special cases treated in lfT9ll20l ). 

So far we have only considered the Euclidean Voronoi and Delaunay diagrams, but a con¬ 
siderable amount of literature exists on Voronoi and Delaunay diagrams under other norms and 
so-called convex distance functions; see Section[2]for details. 

Chew IfTTl showed that the Delaunay triangulation of P under the L\- or Loo-metric experiences 
only a near-quadratic number of events, if the motion of the points of P is algebraic of bounded 
degree. In the companion paper 0, we present a kinetic data structure for maintaining the Voronoi 
diagram and Delaunay triangulation of P under a polygonal convex distance function for an 
arbitrary convex polygon Q. (see Section[2]for the definition) that processes only a near-quadratic 
number of events, and can be updated in O (1 og n) time at each event. Since a regular convex 
/c-gon approximates a circular disk, it is tempting to maintain the Delaunay triangulation under a 
polygonal convex distance function as a (hopefully substantial) portion of the Euclidean Delaunay 
graph of P. Unfortunately, the former is not necessarily a subgraph of the latter |8j|. 

1 We assume the motion of the points to be sufficiently generic, so that no more than four points can become cocircular 
at any given time, and so that equality in (HT is not a local maximum of the left-hand side. 
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Many subgraphs of DT(P), such as the Euclidean minimum spanning tree (MST), Gabriel 
graph, relative neighborhood graph, and a-shapes, have been used extensively in a wide range of 
applications (see e.g. |8]). However no sub-cubic bound is known on the number of discrete changes 
in their structures under an algebraic motion of the points of P of bounded degree. Furthermore, 
no efficient kinetic data structures are known for maintaining them, for unlike DT(P), they may 
undergo a "non-local" change at a critical event; see HUS [ED f° r some partial results on maintaining 
an MST. 

Our results. Stable Delaunay edges: We introduce the notion of a-stable Delaunay edges, for a 
fixed parameter a > 0, defined as follows. Let pq be a Delaunay edge under the Euclidean norm, 
and let A pqr + and A pqr~ be the two Delaunay triangles incident to pq. Then pq is called a-stable 
if its opposite angles in these triangles satisfy 

Zpr + q + Zpr~q < 7 z — a. (2) 

As above, the case where pq lies on conv(P) is treated as if r~, say, lies at infinity, so that the 
corresponding angle Zpr~q is equal to 0. An equivalent and more useful definition, in terms of the 
dual Voronoi diagram, is that pq is a-stable if the angles at which p and q see their common Voronoi edge 
e vq are at least a each. See Figure[2ja). 

In the case where pq lies on conv(P) the corresponding dual Voronoi edge e P q is an infinite ray 
emanating from some Voronoi vertex x. We define the angle in which a point p see such a Voronoi 
ray to be the angle between the segment px and an infinite ray parallel to e p q emanating from p. 
With this definition of the angle in which a point of conv(P) sees a Voronoi ray, it is easy to check 
that the alternative definition of a-stability is equivalent to the ordinal definition also when pq lies 
on conv(P). We call the Voronoi edges corresponding to a-stable Delaunay edges a-long (and call 
the remaining edges a-short). See Figure [2] Note that for a = 0, when no four points are cocircular, 
© coincides with ([TJ. 




Figure 2. (a) An a-stable edge pq in DT(P), and its dual edge e p q = ab in VD(P); the fact that Zr+ + Zr~ = 7r — a. 
follows by elementary geometric considerations, (b) A SDG of the Delaunay triangulation DT(P) shown in Figure[l](a), 
with a = n/8. 

A justification for calling such edges stable lies in the following observation: If a Delaunay edge 
pq is a-stable then it remains in DT(P) during any continuous motion of the points of P for which 
every angle Zprq, for r E P \ {p,q}, changes by at most a/2. This is clear because, as is easily 
verified, at any time when pq is a-stable we have Zpr + q + Zpr q < n — a for any pair of points r + . 
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r~ lying on opposite sides of the line £ supporting pq, so, if each of these angles changes by at most 
a/2 we still have Zpr + q + Zpr~q < n for every such pair r + , r~, implying that pq remains an edge 
of DT(P)J^] Hence, as long as the "small angle change" condition holds, stable Delaunay edges 
remain a "long time" in the triangulation. Informally speaking, the non-stable edges pq of DT(P) 
are those for which p and q are almost cocircular with their two common Delaunay neighbors r + , 
r~, and hence pq is more likely to get flipped "soon." 

Stable Delaunay graph: For two parameters 0 < a < a.' < tt, we call a subgraph 9 of DT(P) 
an (a', a)-stable Delaunay graph (an (a', a)-SDG for short) if 

(51) every edge of 9 is a-stable, and 

(52) every a'-stable edge of DT(P) belongs to 9- 

Note that an (a', a)-SDG is not uniquely defined even for fixed a, a' because the edges that are 
a-stable but not a'-stable may or may not be in 9- Throughout this paper, a' will be some fixed (and 
reasonably small) multiple of a. 



Figure 3. A convex distance function, induced by Q, that is a-close to the Euclidean norm. 


Our main result is that a stable edge of the Euclidean Delaunay triangulation appears as stable 
edge in the Delaunay triangulation under any convex distance function that is sufficiently close 
to the Euclidean norm (see Section [2] for more details). More precisely, we say that the distance 
function induced by a compact convex set Q is a-close to the Euclidean norm if Q is contained in 
the unit disk Do and contains the disk Dj = (cos a)Do both centered in the origin^] See Figure 
[3j In particular, for k = n/a, the regular k- gon Q/ c is such a set, as easy trigonometry shows. We 
prove the following: 

Theorem 1.1. Let P be a set ofn points in R 2 3 , a > 0 a parameter, and Q a compact, convex set inducing a 
convex distance function dg(-, •) that is a-close to the Euclidean norm. Then the following properties hold. 

(i) Every lla-stable Delaunay edge under the Euclidean norm is an a-stable Delaunay edge under dQ. 

(ii) Symmetrically, every lla-stable Delaunay edge under dg is also an a-stable Delaunay edge under the 
Euclidean norm. 

2 This argument also covers the cases when a point r crosses l from side to side: Since each point, on either side of t, 
sees pq at an angle of < n — a, it follows that no point can cross pq itself - the angle has to increase from n — a. to n. Any 
other crossing of £ by a point r causes /Lprq to decrease to 0, and even if it increases to a/2 on the other side of t, pq is 
still an edge of DT, as is easily checked. 

3 The Hausdorff distance between Q and Dq is at most 1 — cos « fs r/2. 
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In particular, if Q is a regular k- gon for k > 'In/a., then the above theorem holds for Q. In the 
companion paper [[3j, we have presented an efficient kinetic data structure for maintaining the 
Delaunay triangulation and Voronoi diagram of P under a polygonal convex distance function. 
Using this result, we obtain the second main result of the paper: 

Theorem 1.2. Let P be a set of n moving points in IR 2 under algebraic motion of bounded degree, and 
let a. G (0,7r) be a parameter. A Euclidean (8a, a)-stable Delaunay graph of P can be maintained by a 
linear-size KDS that processes O ( \nA r (n)) events and updates the SDG at each event in 0(log n) time. 
Here r is a constant that depends on the degree of the motion of the points ofP, and A r (n) is the maximum 
length of a Davenport-Schnizel sequence of order r. 



2 Preliminaries 

This section introduces a few notations, definitions, and known results that we will need in the 
paper. 

We represent a direction in R 2 as a point on the unit circle S 1 . For a direction u G S 1 and an 
angle 6 G [0,27t), we use u + 9 to denote the direction obtained after rotating the vector ou by 
angle 6 in clockwise direction. For a point i GR 2 and a direction u G S 1 , let u[x] denote the ray 
emanating from x in direction u. 

Q-distance function. Let Q be a compact, convex set with non-empty interior and with the origin, 
denoted by o, lying in its interior. A homothetic copy Q' of Q can be represented by a pair (p, A), 
with the interpretation Q' = p + AQ; p is the placement (location) of the center o of Q', and A is its 
scaling factor (about its center). Q defines a distance function (also called the gauge of Q) 

dQ{x,y ) = min{A | y G x + AQ}. 

Note that, unless Q is centrally symmetric with respect to the origin, dq is not symmetric. 

Given a finite point set P C R 2 and a point p G P, we denote by Vor^(p), VD^(P), and 
DT^(P) the Voronoi cell of p, the Voronoi diagram of P, and the Delaunay triangulation of P, 
respectively, under the distance function dg(-, •); see Figure [l] (b). To be precise (because of the 
potential asymmetry of cIq), we define 

Vor(p) = {x G R 2 | dQ(x,p ) < dQ(x,p')\/p' G P}, 

and then VD^(P) and DT®(P) are defined in complete analogy to the Euclidean case. We refer the 
reader to for formal definitions and details of these structures. Throughout this paper, we will 
drop the superscript Q from Vor®, VD®, DT® when referring to them under the Euclidean norm. 

For a point z G R 2 , let Q[z\ denote the homothetic copy of Q centered at z such that its 
boundary touches the Q-nearest neighbor(s) of z in P, i.e., Q[z] is represented by the pair (z. A) 
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Figure 4. (a) Homothetic copies Q p («), centered at c € u[p] (dashed), and Q pq (u) (solid) of Q. The ray u[p] hits bfy at 
the center c pq of Q pq {u), which is defined, when Q is smooth at its contact with p, if and only if q and Q p (u ) lie on the 
same side of £ p (u); (b) Sliding Q pr (u) away from q. 


where A = min p6 p <1q ( z,p ). In other words, Q [z] is the largest homothetic copy of Q that is centered 
at z whose interior is P-empty. We also use the notation Q p (u ) to denote a "generic" homothetic 
copy of Q which touches p and is centered at some point on u[p]. See Figure [1] (a). Note that all 
homothetic copies of Q p {u) touch p at the same point £ of 9 Q, and therefore share the same tangent 
at p. This tangent is unique if Q is smooth at £, and we denote it by £ p (u). If Q is not smooth at £ 
then there is a nontrivial range of tangents (i.e., supporting lines) to Q at l,; we can take £ p (u) to be 
any of them, and it will be a supporting line to all the copies Q p (u ) of Q. 

For a pair of points p,q E P, let bfy denote the Q-bisector of p and q —the locus of all placements 
of the center of any homothetic copy Q' of Q that touches p and q. If Q is strictly convex or if Q 
is not strictly convex but no two points are collinear with a straight segment on 9 Q, then b ( pq is a 
one-dimensional curve and any ray u [p] that hits b® q does so at a unique point. For such a direction 
u and a pair of points p,q E P, let Q pq (u) denote the homothetic copy of Q that touches p and q, 
whose center is u[p] n b® q . 

If Q is not strictly convex and p, q are points in P such that pq is parallel to a straight portion e 
of 9 Q then b pq is not one-dimensional. In this case Q pq (u ) is not well defined when u is a direction 
that connects e to the center c of Q. As is easy to check, in any other case the ray u[p] either hits b pq 
at a unique point which determines Q pq {u), or entirely misses bj? q . See the companion paper (3) for 
a detailed discussion of this phenomenon. 

A useful property of the Q-bisectors is that any two bisectors b® q , b ( j; r with a common generating 
point p, intersect exactly once, namely, at the center c* of the unique homothetic copy of Q that 
simultaneously touches p, q and r Ifl7l . For this property to hold, though, we need to assume that (i) 
the points p, q, r are not collinear, and (ii) either Q is strictly convex, or, otherwise, that none of the 
directions pq, pr is parallel to a straight portion of 9 Q. (The precise condition is that b pq and b pr be 
one-dimensional in a neighborhood of c*.) The local topology of the restricted Q-Voronoi diagram 
VD < '-(\p,q, r}) near c* is largely determined by the orientation of the triangle A pqr. Specifically, 
assume with no loss of generality that pr is counterclockwise to pq, and let u* be the direction of the 
ray pc*. Refer to Figure Mhb). If we continuously rotate a ray u[p], for u E S 1 , in counterclockwise 
direction from u* [p], the corresponding copy Q pr (u) will slide away from its contact with q because 
the portion of Q pr (u ) to the right of pr shrinks during the rotation. Therefore, the rotating ray u[p\ 
either misses b entirely or hits bfy after b ( pr . A symmetric phenomenon, with q and r interchanged, 
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takes place if we rotate the ray u[p] in clockwise direction from u* [p\. 

It is known that Vor®(p), for every point p E P, is star-shaped (8), which implies that each 
Voronoi edge is fully enclosed between the two rays that emanate from p, or from q, through its 
endpoints. We remark that, unlike the Euclidean case, the angles Z xpxj , Zxqy need not be equal in 
general. 

Finally, we extend the notion of stable edges to DT^(P). We call an edge pi] E DT® (P) a.-stable 
if the following property holds for the dual edge in VD^(P): 

Each of the points p, q sees their common Q-Voronoi edge eat angle at least a. That is, if x 
and y are the endpoints ofethen min{Zxpy, Zxqyj > a. 

This definition coincides with the definition of a-stability under the Euclidean norm when Q is the 
unit disk (and in this case both angles are equal). 


Remark. If Q is not strictly convex (and pq is parallel to a straight portion of 9Q), the endpoints 
~pq may be not well defined. In this case, we resort to the following, more careful definition of 


of e9 


a-stability. A ray u[p\ is said to (properly) cross efy only if the copy Q P q(u) is uniquely defined. 
The center of such a copy Q pq {u) necessarily lies within the one-dimensional portion e ( p[] of Q pq (u), 


which is easily seen to be non-empty and connected. We say that the edge is a-stable if the set of 
rays u[p] properly crossing e pq spans an angle of at least a, and a symmetric condition holds for 
the rays emanating from q. In other words, our notion of ^-stability ignores the two-dimensional 
regions of (if these exist) 


Polygonal convex distance function. As mentioned in the introduction, we will be considering 
the case when Q is a regular k- gon, for some even integer k > 2n/a., centered at the origin. Let 
Vo, ■ ■ ■, Ujt-i be its sequence of vertices arranged in clockwise direction. For each 0 < j < k, let 
Uj be the direction of the vector that connects Vj to the center of Q (see Figure [5] (a)). We will use 
bpq, Vor°(p),VD°(P),DT <> (P) to denote b p q ,Vor®(p),VD^(P),DT®(P), respectively, when Q is a 
regular k-gon. 

We say that P is in general position (with respect to Q) if no three points of P lie on a line, no 
two points of P lie on a line parallel to an edge or a diagonal of Q, and no four points of P are 
Q-cocircular , i.e., no four points of P lie on the boundary of a common homothetic copy of Q. 

The placements on b pq at which (at least) one of p and q, say, p, touches Q' at a vertex is called a 
corner placement (or a corner contact) at p; see Figure |5](b). We also refer to these points on b pp as 
breakpoints. We call a homothetic copy of Q whose vertex Vj touches a point p, a vj-placement ofQ at 
V- 

The following property of b pcj is proved in [[3} Lemma 2.5]: 


Lemma 2.1. Let Q be a regular k-gon, and let p and q be two points in general position with respect to Q. 
Then b pc) is a polygonal chain with k — 2 breakpoints and the breakpoints along b pq alternate between corner 
contacts at p and corner contacts at q. 

4 As is easy to check, the one-dimensional portion e p q of efy varies continuously (in Hausdorff sense) with any 
sufficiently small perturbation of p and q within P. Furthermore, it is the only such portion: If a ray u [p] hits e p q outside 
(i.e., within its two-dimensional portion), there is a symbolic perturbation of p and q causing u[p] to completely miss 
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(a) 


(b) 


Figure 5. (a) A regular octagon Q centered at the origin, (b) The bisector b pq for the regular octagon Q; it has six 
breakpoints and the corner contacts along b pq alternate between contacts at p (hollow circles) and contacts at q (filled 
circles). 


3 DT°(P) and Euclidean SDG's 


In this section we first prove a slightly stronger version of Theorem 1.1 for the case when Q is a 
regular fc-gon for some even integer k > 2n/a, and then prove Theorem 1.2 We follow the notation 
in Section]!} and, for simplicity, we assume that a = 2n/k. 



Figure 6. (a) The function cpj [p, q\, which is equal to the radius of the circle that pass through p and q and whose center 
lies on Uj [p]. (b) The bisector b pcf and the function cpj [p, q], which is equal to the radius of the circle that circumscribes 
the py-placement of Q at p that also touches q. (c) The case when cpj[p, q] = oo while cpj[p, q] < oo. In this case q lies in 
one of the shaded wedges. 


For any pair p,q E P, let cpj [p, q] denote the distance from p to the point Uj [p] (1 b pq ) we put 
< pj[p,q] = oo if itjlp] does not intersect b pq . See Figure[6](a). The point q minimizes (pi[p,q'\, among 
all points q ' for which Uj [p] intersects b pq t, if and only if the intersection between b pq and w, [p\ lies 
on the Voronoi edge e pq . We call q the neighbor of p in direction n lr and denote it by /V, (p). 

Similarly, let cpj[p , q\ denote the distance from p to the point Uj [p] PI b pcj ; we put cp - [p, q\ = oo if 
Uj [p] does not intersect b pq . If cp- [p, q] < oo then the point b pcj n Uj [p] is the center of the uy-placement 
Q' of Q at p that also touches q, and there is a unique such point. The value cpj [p, q] is equal to the 
circumradius of Q'. See Figure |6j(b). The neighbor NJ [p\ of p in direction uj is defined to be the 
point q G P \ {p} that minimizes cpj[p,q]- 

Note that (pj[p,q\ < oo if and only if the angle between pq and Uj[p] is smaller than n/2. In 
contrast, cpj [p, q] < oo if and only if the angle between pq and Uj[p\ is at most n/2 — n/k = 
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7 t /2 — a/2. Moreover, we have <Pj[p,q\ < [p, q] (see Figure |6l). Therefore, < °o always 

implies <p ; [p, q] < oo, but not vice versa; see Figure[6](c). Note also that in both the Euclidean and 
the polygonal cases, the respective quantities Nj [p] and Nj [p\ maybe undefined. 


Lemma 3.1. Let p,q G P be a pair of points such that Nj(p) = q for h > 3 consecutive indices, say 
0 < j < h — 1. Then for each of these indices, except possibly for the first and the last one, we also have 


Nf[p] = q. 


Proof. Let w\ (resp., wf) be the point at which the ray uo[p\ (resp., u h _-\ [p\) hits the edge e pq in 
VD(P). (By assumption, both points exist.) Let D\ and D 2 be the disks centered at W\ and 102, 
respectively, and touching p and q. By definition, neither of these disks contains a point of P in 
its interior. The angle between the tangents to D 1 and D 2 at p or at q (these angles are equal) is 
f> = (h — l)a; see Ligure[7](a). 





Figure 7. (a) The angle between the tangents to Dj and D 2 at p (or at q) is equal to = jB = (h — l)a. (b) The 

disks D and D+ and the homothetic copy Qj of Q; t! n D = qq' C e’. (c) £' forms an angle of at least a. /2 with each of the 
tangents to Dj, Dt at q, and the edge e’ = C Dj U D 2 . 


Lix an arbitrary index 1 < / < h — 2, so Uj[p] intersects e p q and forms an angle of at least a 
with each of pw v pw 2 . Let Q ( := Qpq(uj ) be the ^-placement of Q at p that touches q. To see that 
such a placement exists, we note that, by the preceding remark, it suffices to show that the angle 
between pq and Uj[p\ is at most tc /2 — a/2; that is, to rule out the case where q lies in one of the 
shaded wedges in Ligure[6](c). This case is indeed impossible, because then one of Uj-i[p], U i+M 
would form an angle greater than 7r/2 with pq, contradicting the assumption that both of these 
rays intersect the (Euclidean) b p q. The claim now follows from the next lemma, which shows that 
QjC. Di U D 2 , which implies that int Qj FI P = 0 and thus NJ [p] = q, as claimed. 'jpj| 


Lemma 3.2. In the notation in the proof of Lemma 3.1 Qj C Di U D 2 ,for 1 < / < h — 2. 


Proof. Lix a value of l<j<h-2. Let e' be the edge of Qj passing through q; see Ligure[7](b). 
Let D be the disk whose center lies on Uj[p] and which passes through p and q, and let D + be the 
circumscribing disk of Qj. Since p G 9D IT 9D + , q G 9 D IT intD + , and D and D + are centered on 
the ray Uj [q] emanating from p, it follows that D C D + . The line f containing e' crosses D in a 
chord qq’ that is fully contained in e’, as qq’ = D IT £' C D + IT f = e’. 

The angle between the tangent to D at q, denoted by r, and the chord qq’ is equal to the angle at 
which p sees qq'. This angle is smaller than the angle at which p sees e', which in turn is equal to 
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a/2. Recall that u q [p] makes an angle of at least a with each of pw-\ and piv 2 , therefore r forms an 
angle of at least a with each of the tangents to D \, D 2 at q. Combining this with the fact that the 
angle between r and e' is at most a/2, we conclude that e' forms an angle of at least a/2 with each 
of these tangents; see Figure [7] (c). 

The line (! marks two chords q\q,qq 2 within the respective disks Dj, D 2 . We claim that e' is 
fully contained in their union i] \ q 2 - Indeed, the angle i] \ pt] is equal to the angle between P and the 
tangent to D 1 at q, so Gq\ pt] > a/2. On the other hand, the angle at which p sees e' is a/2, which is 
no larger. This, and the symmetric argument involving D 2 , are easily seen to imply the claim. 

Now consider the circumscribing disk D + of Q ( . Denote the endpoints of e' as ci \ and a 2 , where 
fli lies in i]]i] and a 2 lies in qq 2 - Since the ray pa ] hits 9D + before hitting 9D 1 , and the ray pq hits 
these circles in the reverse order, it follows that the second intersection of 9Di and 9D + (other than 
p) must lie on a ray from p which lies between the rays pa ,, pq and thus crosses e'. See Figure [71(c). 
Symmetrically, the second intersection point of 9D2 and 9D + also lies on a ray which crosses e. It 
follows that the arc of 9D + delimited by these intersections and containing p is fully contained in 
D 1 U D 2 . Flence all the vertices of Q ; (which lie on this arc) lie in Di U D 2 . This, combined with the 
fact, established in the preceding paragraph, that e' C q 1 q 2 implies that Q ; C Di U D 2 . Jpl 


Next, we use Lemma 3.1 to prove its converse. Specifically, we prove the following lemma. 


Lemma 3.3. Let p,q G P be a pair of points such that N?[p] = q for at least five consecutive indices 
j G {0,... ,k — 1}. Then for each of these indices, except possibly for the first tzvo and the last two indices, 
we have Nj [p] = q. 


Proof Again, assume with no loss of generality that N'j’lp] = q for 0 < j < h — 1, with h > 5. 
Suppose to the contrary that, for some 2 < j < h — 3, we have Nj[p] fi q. By assumption, Nf [p] = q, 
for each 0 < i < h — 1, and therefore we have cf>i[p, q\ < (p- [p, q] < 00 , for each of these indices. 
In particular, we have cpj[p, q] < [p, q\ < 00 , so there exists r G P for which cpj[p, r\ < <pj[p,q\. 

Assume with no loss of generality that r lies to the left of the line from p to q. In this case we claim 
that <P/-i [p,r] < cpj-i [p,q] < 00 and <p ,_2 [p, r] < (pj- 2 [p,q\ < 00 . 

Indeed, the boundedness of (pj 1 [p, q\ and cpj _ 2 [p, q] has already been noted. Moreover, because 
r lies to the left of the line from p to q, the orientation of b pr lies counterclockwise to that of b pq . 
This, and our assumption that ufip] hits b pr before hitting b pq , implies that the point b pr fl b pq 
lies to the right of the (oriented) line through Uj[p]; see Figure |8l Flence, any ray p emanating 
from p counterclockwise to Uj[p\ that intersects b pq must also hit o pr before hitting b pq , so we have 
(Pj-l[p,r\ < cpj-i [p, q\ and cpj_ 2 [p,r] < cpj- 2 [p,q] (since j > 2, both w ; _i[p] and Uj- 2 [p] intersect b pq ), 
as claimed. Now applying Lemma |3d| to the point set {p, q,r} and to the index set { / — 2 ,j — 1 ,j}, 
we get that (pfi [p, r] < q)J _ } [p, q]. But this contradicts the fact that NJ_ q [p] = q. The case where r 
lies to the right of pq is handled in a fully symmetric manner, using the indices { /, / + l, j + 2}. Q 

Combining Lemmas |3.1| and |3.3[ we obtain the following stronger version of Theorem |1.1| 

Theorem 3.4. Let P be a set ofn points in R 2 , a > 0 a parameter, and Q a regular k-gon with k > 2n/a.. 
Then the follozving properties hold. 

(i) Every 4 oc-stable edge in DT(P) is an a-stable edge in DT <> (P). 

(ii) Every 6a-stable edge in DT*(P) is also an oc-stable edge in DT(P). 
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Figure 8. Proof of Lemma 3.3 


Proof. Let pq be a 4a-stable edge of DT(P). Then the corresponding edge e pq in VD(P) stabs at least 
four rays Uj [p] emanating from p, and, by Lemma 3.1 Nj [p] = q for at least two of these values of j. 
Therefore, p sees the edge e pq 


,0 in VD*(P) at an angle at least a. Similarly q sees the edge et„ at an 


angle at least a. Conversely, if pq is 6a-stable in DT°(P) then e pq meets at least six rays Uj[p], and 
then Lemma 3.3 is easily seen to imply that p (and, symmetrically, q too) sees e pq at an angle at least 

DC. □ 


The next lemma gives a slightly different characterization of stable edges, which is more 
algorithmic and which will be useful in maintaining a SDG under a constant-degree algebraic 
motion of the points of P. 

Lemma 3.5. Let 9 be the subgraph o/DT°(P) composed of the edges whose dual Q-Voronoi edges contain 
at least eleven breakpoints. Then 9 is an (8a, a)-SDG ofP (in the Euclidean norm). 


Proof. Let p, q G P be two points. If (p, q) is an 8a-stable edge in DT(P) then the dual Voronoi edge 
e pq stabs at least eight rays Uj[p] emanating from p, and at least eight rays Uj[q\ emanating from q. 
Lemma |3 .1 1 implies that VD°(P) contains the edge e pq with at least six breakpoints corresponding 
to corner placements of Q at p that touch q, and at least six breakpoints corresponding to corner 
placements of Q at q that touch p. Therefore, e pq contains at least twelve breakpoints, so ( p,q ) G 9- 
Conversely, suppose p,q G P define an edge e pq in VD <> (P) with at least eleven breakpoints. By 
the interleaving property of breakpoints, stated in Lemma |2.1[ we may assume, without loss of 
generality, that at least six of these breakpoints correspond to empty corner placements of Q at p 
that touch q. Lemma 3.3 implies that VD(P) contains the edge e pq , and that this edge is hit by at 
least two consecutive rays Uj[p]. But then the ( p,q ) is a-stable in DT(P). □ 


In a companion paper [j3|, we describe a kinetic data structure (KDS) for maintaining DT <> (P). 
As shown in that paper, it can also keep track of the number of breakpoints for each edge of DT <> (P). 
If each point of P moves along an algebraic trajectory of bounded degree, then the KDS processes 
0(k 4 n\ r (n)) events, where r is a constant depending on the complexity of motion of P. A change 
in the number of breakpoints in a Voronoi edge is an event that the KDS can detect and process. 
As discussed in detail in (3|, many events, so-called singular events, that occur when an edge of 
DT <> (P) becomes parallel to an edge of Q, can occur simultaneously. Nevertheless, each of the 
events can be processed in 0(log n ) time, and their overall number is within the bound cited above. 
We maintain the subgraph 9 of DT°(P), consisting of the edges of DT <> (P) that have at least eleven 
breakpoints, which, by Lemma [3.5[ is an (8a, a)-Euclidean SDG. Putting everything together, we 
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obtain a KDS that maintains an (8a, a)-SDG of P. It uses linear storage, it processes O (-^}iA r (n) j 
events, where r is a constant that depends on the degree of the motion of the points of P, and it 
updates the SDG at each event in 0(log zz) time. This proves Theorem 

Remarks, (i) We remark that DT°(P) can undergo O(n) changes at a time instance when Q(zz) 
singular events occur simultaneously, say, when pq becomes parallel to an edge of Q, but all these 
changes occur at the edges incident to p or q in DT(P). However, only 0(1 /a) edges among 
them can have at least eleven breakpoints, before or after the event. Hence, 0(1/a) edges can 
simultaneously enter or leave the (8a, a)-SDG S of Theorem |1.2| 

(ii) Note that there is a slight discrepancy between the value of k that we use in this section 
(k > In /a), and the value needed to ensure that the regular /r-gon Qp is a-close to the Euclidean 
disk, which is k > n/ a. This is made for the convenience of presentation. 

(iii) An interesting open problem is whether the dependence on a can be improved in the above 
KDS. We have developed an alternative scheme for maintaining stable (Euclidean) Delaunay 
graphs, which is reminiscent of the kinetic schemes used by Agarwal et al. 0) for maintaining 
closest pairs and nearest neighbors. It extracts a nearly linear number of pairs of points of P 
that are candidates for being stable Delaunay edges and then sifts the stable edges from these 
candidate pairs using the so-called kinetic tournaments ItTOll . Although the overall structure is 
not complicated, the analysis is rather technical and lengthy, so we omit this KDS from this 
version of the paper; it can be found in the arXiv version 0. In summary, the resulting KDS is 
of size 0((n/a 2 ) logzz), it processes a total of 0((n/a.) 2 (5 r {n/ a) log 2 zzlog(log(zz)/a)) events, and 
it takes a total of 0((zz/a) 2 log 2 zt(log 2 zz + f5 r (n/cc) log 2 zzlog 2 (log(zz)/a))) time to process them; 
here j8,-(zz) = A r (zz)/zz is an extremely slowly growing function for any fixed r. The worst-case time 
of processing an event is O(log 4 (zz) /a)). Another advantage of this data structure is that, unlike the 
above KDS, it is local in the terminology of flOl . Specifically, each point of P is stored, at any given 
time, at only 0(log 2 (zz) /a 2 ) places in the KDS. Therefore the KDS can efficiently accommodate an 
update in the trajectory of a point. 

4 Stability under Nearly Euclidean Distance Functions 

In this section we prove Theorem [hT] for an arbitrary convex distance function that is a-close to the 
Euclidean norm (see the Introduction for the definition). 

Let Q be a compact, convex set that contains the origin in its interior, and let dq denote the 
distance function induced by Q. Assume that Q is a-close to the Euclidean norm. 

For sake of brevity, we carry out the proof assuming that Q is strictly convex, i.e., the relative 
interior of the chord connecting two point x and y on 9Q is strictly contained in the interior of Q 
(there are no straight segments on the boundary of Q). The proof also holds verbatim when Q is 
not strictly convex, provided that no pair of points p, q E P is such that pq is parallel to a straight 
portion of 9 Q. 

We assume that P is in general position with respect to Q, in the sense that no three points of P 
lie on a line, and no four points of P are Q-cocircular, i.e., no four points of P lie on the boundary of 
a common homothetic copy of Q. 

Recall that for a direction u and for a point p E P, Qp(u) denotes a "generic" homothetic copy of 
Q that touches p and is centered at some point on u[p\. See Figure |I](a). As mentioned in Section[2j 
all homothetic copies of Qp(u ) touch p at (points corresponding to) the same point £ E 9 Q and 


1.2 


12 



therefore share the same tangent £ v {u) at p. If Q is not smooth at £, there is a range of possible 
orientations of such tangents. In this situation we let £p(u) denote an arbitrary tangent of this kind. 
In addition, if Q is smooth at £ then Q P q(u) exists for any point q that lies in the same side of£ p (u) 
as Qp(u). Otherwise, this has to hold for every possible tangent £p(u) at £, which is equivalent to 
requiring that q lies in the wedge formed by the intersection of the two halfplanes bounded by the 
two extreme tangents at £ and containing Q v (u). 

Remarks. (1) An important observation is that, when q satisfies these conditions, Qpq(u) is unique, 
unless all the three following conditions hold: (i) Q is not strictly convex, (ii) pq is parallel to 
straight portion e of dQ, and (iii) u is a direction connecting some point on e to the center c of Q. 
We leave the straightforward proof of this property to the reader. The proof of the theorem exploits 
the uniqueness of Q P q(u), and breaks down when it is not unique. In fact, this is the only way in 
which the assumptions concerning strict convexity are used in the proof. 

(2) With some care, our analysis applies also if (the directions of) some pairs pq are parallel to 
straight portions of dQ, in which case Q p q(u) is not uniquely defined for certain directions it. 
This extension requires the more elaborate notion of a-stability, which ignores the possible two- 
dimensional portions of ejy see Section 2 for more details. Informally, this allows us to avoid the 
"problematic" directions u in which Qpq[u) is not unique. (The latter happens exactly when u[p\ 
hits \>pq within one of its two-dimensional portions.) We note, though, that the loss in the amount 
of stability caused by ignoring a two-dimensional portion of is at most 2 a. which is an upper 
bound on the angular span of directions that connect a straight portion of dQ to its center. This 
latter property holds since Q is a-close to the Euclidean disk; see below for more details. 

The proof of Theorem 1 1.1 1 relies on the following three simple geometric properties. Recall that 
the ^-closeness of Q to the Euclidean norm means that Dj = (cos a ) Dq CQC Dq. 




Figure 9. Illustrations for: (a) Claim |4~l) (b) Claim |4~2| 


Claim 4.1. Let x be a point on dQ, and let £ be a supporting line to Q at x. Let y be the point on 9Dq 
closest to x (x and y lie on the same radius from the center o), and let 7 be the arc ofdDp that contains y 
and is bounded by the intersection points of £ with d Dq. Then the angle between t and the tangent, r, to Dq 
at any point along 7 , is at most a. 

Proof. Denote this angle by 6. Clearly 9 is maximized when r is tangent to Dq at an intersection of 
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£ and dD 0 ; see Figure [9] (a). For this value of 0, it is easy to verify that the distance from o to £ is 
cos 9. But this distance has to be at least cos a, because £ fully contains Dj = (cos a)Do on one side. 
Hence cos 9 > cos a, and thus 9 < a, as claimed. □ 


Remark. An easy consequence of this claim is that the angle in which the center of Q sees any 
straight portion of dQ (when Q is not strictly convex) is at most 2a. 


Claim 4.2. Let x and y be two points on dQ, and let L x and £ y be supporting lines of Q at x and y, 
respectively. Then the difference between the (acute) angles that £ x and £ y form zvith xy is at most 2a. 


Proof. Denote the two angles in the claim by 9 X and 9 y , respectively. Continue the segment xy 
beyond x and beyond y until it intersects Do at x' and i/, respectively. Let T\ and T 2 denote the 
respective tangents to D () at x' and at y'. See Figure [9] (b). Clearly, the respective angles 9\, 02 


4.1 


between the chord x'y' of Do and r\, T 2 are equal. By Claim 
to T 2 and £ y ) we get that 1 9\ — 9 X | ^ a and |02 — 9 y \ < a, and the claim follows. 


(applied once to X\ and i x and once 

□ 


Do 



Figure 10. Proof of Claim |4~3| 


Claim 4.3. For a point p G dQ, any tangent £ p to Q at p forms an angle at most a zvith any line orthogonal 
to op. 


Proof. See Figure [lOj Consider the chord £ = £ p fl Dq, and let 7 denote the arc of dD 0 determined by 
£ and containing the intersection point z of dD ( j and the ray op. By Claim 4.1 the angle between £ p 
and the tangent to Dq at z is at most a. Since this tangent is orthogonal to op, the claim follows. □ 


Remark. Clearly, Claims 4.1 -4.3 continue to hold for any homothetic copy of Q, with a correspond¬ 
ing translation and scaling of Dq and D 7 . 


Let Q pq (u ) (resp., Q pcj (u)) denote the portion of Q P q(u) that lies to the left (resp., right) of the 
directed line from p to q. Let D p q(u) denote the disk that touches p and q, and whose center lies on 
u[p]. 

We next establish the following lemma, whose setup is illustrated in Figure [Tl] (a). It provides 
the main geometric ingredient for the proof of Theorem |1.1| 

Lemma 4.4. (i) Let u G S 1 be a direction such that both Q py (u) and D pq (u + 5a) are defined. Then the 
region Q^ q (u) is fully contained in the disk D pq (u + 5a). 

(ii) Let u G S 1 be a direction such that both Q pq {u) and D pq (u — 5a) are defined. Then the region Q pq (u) 
is fully contained in the disk D pq (u — 5a). 
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Proof. It suffices to establish Part (i) of the lemma; the proof of the other part is fully symmetric. 

Refer to Figure [TT| (b). Let £ p = £ p {u) be any supporting line of Q P q(u) at p, as defined above, 
and let T p = r p (u) be the line through p that is orthogonal to u (which is also the tangent to 
Dpq(u)). By Claim 4.3 the angle between £ p and r p is at most a. We next consider the tangent T p to 


D pq (u + 5a) at p. Since the angle between r p and is 5a, it follows that the angle between £ p and 
T p is at least 4a. The preceding arguments imply that, when oriented into the right side of pq, £ p 
lies between pq and t+, and the angle between £ p and r+ is at least 4a. 

This implies that, locally near p, 3Q+ ; (m) penetrates into D p[ fu + 5a). This also holds at q. 
To establish the claim for q (which is not symmetric to the claim for p, because the center c of 
Q P q(u) lies on the ray u[p\ emanating from p, and there is no control over the orientation of the 
corresponding ray qc emanating from q), we note that, by Claim |4.2| the angles between pq and 
any pair of tangents £ p , £ q to Q pq (u) at p, q, respectively, differ by at most 2a, whereas the angles 
between pq and the two tangents T p , T+ to D pq (u + 5a) at p, q, respectively, are equal. This, and 
the fact that the angle between and l p is at least 4a, imply that, when oriented into the right side 
of pq, £q lies between cfp and t+, which thus implies the latter claim. Note also that the argument 
just given ensures that the angle between t ' and £ q is at least 2a. 




Figure 11. (a) The setup of Lemma 4.4 (b) Proof of Lemma 4.4 (i): 3Q+ (u) cannot cross dD p q(u + 5a) at any third point 


It therefore suffices to show that 9Q+(u) does not cross dD pq (u + 5a) at any third point (other 
than p and q). Suppose to the contrary that there exists such a third point zu, and consider the 
tangents £ w to Q p q(u) at zv, and rf to D pp (u + 5a) at zv. Consider the two points p and zv, and 
apply to them an argument similar to the one used above for p and q. Specifically, we use the facts 
that (i) the angles between pw and r p , rf are equal, (ii) the angles between pzv and £ p , £ w , for any 
tangent £ w to Q pp (u) at zv, differ by at most 2a, and (iii) the angle between £ p and is at least 4a, 
to conclude that, when oriented into the left side of zvp, £ w lies strictly between zvp and r+. See 
Figure [IT] Similarly, applying the preceding argument to q and w, we now use the facts that (i) the 
angles between zvq and T+, are equal, (ii) the angles between qzv and £ q , £ w differ by at most 2a, 
and (iii) the angle between rf and £ q is at least 2a, to conclude that, when oriented into the right 
side of zvq, £ w lies between zvq and rf or coincides with r+. This impossible configuration shows 
that zv cannot exist, and consequently that Q+ (u) C D pq (u + 5a). □ 
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Proof of Theorem |1.1| — Part (i). Let pq be an 1 la-stable edge in the Euclidean Delaunay triangu¬ 
lation DT(P). That is, the Euclidean Voronoi edge e p q is hit by two rays it [p], u 1 [p] which form an 
angle of at least 11a between them (where u + is assumed to lie counterclockwise to u). Clearly, e p q 
is also hit by any ray u [p] whose direction u belongs to the interval (u , u + ) C S 1 . Let u [p] be such 
a ray whose direction u belongs to the interval (u~ + 5a, u + — 5a) (of span at least a). That is, u[p] 
hits e pq "somewhere in the middle", so all the three disks D pq (u — 5a), D p q(u ) and D pq (u + 5a) are 
defined and contain no points of P in their respective interiors. (Actually, D pq (u) is contained in 
D p q{u — 5a) U D pq {u + 5a), as is easily checked.) 

We next consider the Q-Voronoi diagram VD® (P). We claim that the corresponding edge e p q 
exists and is also hit by u[p]. Since this holds for every u E (u + 5a, u + — 5a), it follows that (p, q) 
is a-stable in VD®(P). 

To establish this claim, we prove the following two properties. 

(i) the homothetic copy Q pq (u) exists, and 

(ii) it contains no points of P in its interior. 


Proof of (i): Assume to the contrary that the copy Q pq {u) is undefined. Consider the respective 
tangents £ p (u) and T p (u) to Q p {u) and D pq (u) at p, where £ p (u) is any tangent to Q p (u) at p that 
separates q from Q p (u)' r such a tangent exists if and only if Q pq (ii) is undefined. (As noted before, 
£p(u) does not depend on the location of the center c of Q on u[p].) By Claim |4~3j the angle between 
£p(u) and r p (u) is at most a. Since Q pq (u) is undefined, the choice of £ p (u) guarantees that q lies 
inside the open halfplane h p {u) bounded by t p (u) and disjoint from u[p]. 

Let p p {u + 5a) (resp., p p (u — 5a)) denote the open halfplane bounded by T p (u + 5a) (resp., 
r p ( u — 5a)) and disjoint from the disk D pq (u + 5a) (resp., D pq (u — 5a)). Since each of the lines 
T p ( u + 5a), T p (u — 5a) makes an angle of at least 5a with t p (m), the halfplane h p (u) supported 
by £p(u) is contained in the union }i p (u + 5a) U p p (u — 5a). Since q is contained in h p (u), at least 
one of these latter halfplanes, say p p (u + 5a), must contain q. However, if q E }i p (u + 5a), the 
corresponding copy D pq (u + 5a) is undefined, a contradiction that establishes (i). 


4.4 


i) implies that Q pq (u) C 
(u) is also 


Proof of (ii): Since both Q pq (u) and D pq (u + 5a) are defined. Lemma 
D pq (u + 5a). Moreover the interior of D pq (u + 5a) is P-empty, so the interior of Q pq 
P-empty. A symmetric argument (using Lemma |L4|ii)) implies that the interior of Q~ q (u) is also 
P-empty. 

This completes the proof of part (i) of Theorem |1.1[ 


Proof of Theorem |l.l| - Part (ii). We fix a direction u E S 1 for which all the three copies Q pq (u), 
Q pq (u — 5a), and Q pq (u + 5a) are defined and have P-empty interiors. Again, Q pq (u) C Q pq (u — 
5a) U Qpq (u + 5a). Since pq is lla-stable under cIq, there is an arc on S 1 of length at least a, so that 
every u in this arc has this property. We need to show that, for every such u, 

(i) the copy D pq (u) is defined, and 

(ii) its interior is P-empty. 

Similar to the proof of Part (i), this would imply that the ray u[p\ hits the edge e pq of VD(P) for 
every u in an arc of length a, so pq is a-stable in DT(P), as claimed. 
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Proof of (i): Assume to the contrary that D pq {u) is undefined, so the angle between the vectors pq 
and u is at least n/1. Let £ p (u — 5a), £p(u), and £ p {u + 5a) be any triple of respective tangents to 
Qpq(u — 5a), Qpq(u), and Q pq {u + 5a) at p. Let h p {u — 5a) (resp., h p {u + 5a)) be the open halfplane 
supported by £ p {u — 5a) (resp., £ p (u + 5a)) and disjoint from Q P q{u — 5a) (resp., Q pq {u + 5a)). 
Claim 4.3 implies that each of the lines £ p {u — 5a), £ p (u + 5a) makes with T p (u), the line orthogonal 
to u[p\ at p, an angle of at least 4a (and at most 6a). Indeed, the claim implies that the angle 
between £ p {u — 5a) and the line T p {u — 5a), which is orthogonal to {u — 5a) [p] at p, is at most a. 
Since the angle between r p (u — 5a) and r p (u) is 5a, the claim for £ p {u — 5a) follows. A symmetric 
argument establishes the claim for £ p (u + 5a). Therefore, the halfplane }tp{u), supported by T p {u) 
and containing q, is covered by the union of h p {u — 5a) and h p {u + 5a). We conclude that at least 
one of these latter halfplanes must contain q. However, this contradicts the assumption that both 
copies Qpq(n — 5a), Q pq {u + 5a) are defined, and (i) follows. 

Proof of (ii): Assume to the contrary that D pCj {u), whose existence has just been established, 
contains some point r of P in its interior. That is, the ray u[p] hits b pr before b p q. In this case D pr {u) 
also exists. With no loss of generality, we assume that r lies to the left of the oriented line from p to 
q. 

We claim that the homothetic copy Q vr {u — 5a) exists and contains q. Indeed, since Q pq {u — 5a) 
exists and is P-empty, it follows that (u — 5a) [p] either hits b ( pr after (in which case the claim 
obviously holds) or misses b ( pr altogether. Suppose that (u — 5a) [p] misses b ( pr . As argued earlier, 
this means that there exists a tangent £ p {u — 5a) to Q p {u — 5a) at p, such that r lies in the open 
halfplane h p {u — 5a) supported by £ p {n — 5a) and disjoint from Q v {u — 5a). 

as before we get that the tangent t p (w) to D p q{u) (at p) to the left of pq is 


By applying Claim 


4.3 


between pq and £ p (u — 5a) and makes with £ p {u — 5a) an angle of at least 4a. It follows that the 
wedge formed by the intersection of £p(u — 5a) and the halfplane to the left of pq is fully contained 
in the halfplane }i p {u ) that is supported by T p (u) and disjoint from D p(J (w); see Figure [l2| (a). But 
then Dp r (u ) is undefined, a contradiction that implies the existence of Q pr {u — 5a). 




(b) 


Figure 12. Proof of Theorem 1.1 ii): (a) Arguing (by contradiction) that Q pr {u — 5 a) exists; (b) the copy Qp r (u — 5a) is 
defined and contains q. Hence, the disk D pr (u), which covers Qf r (u — 5a), must contain q as well. 

We can now assume that Q pr (u — 5a) is defined and contains q. More precisely, q lies in the 
portion Q pr (w — 5a) of Q pr {u — 5a), since q lies to the right of the oriented line from p to r. However, 
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Lemma 


4.4 i), applied to u — 5a, implies that Q pr (u — 5a) is contained in D pr (u), so D pr (u) also 


contains q; see Figure 12] (b). It is however impossible for D pr (u ) to contain q and for D pL] (it) to 
contain r. This contradiction concludes the proof of part (ii) of Theorem |1.1| □ 


5 Properties of stabe Delaunay graphs 

We establish a few useful properties of stable Delaunay graphs in this section. 


Near cocircularities do not show up in an SDG. Consider a critical event during the kinetic 
maintenance of DT(P), in which four points a, b, c, d become cocircular, in this order, along their 
circumcircle, with this circle being P-empty. Just before the critical event, the Delaunay triangulation 
contained two triangles formed by this quadruple, say, abc, acd. The Voronoi edge e ac then shrinks to 
a point (namely, to the circumcenter of abed at the critical event), and, after the critical cocircularity, 
is replaced by the Voronoi edge e^, which expands from the circumcenter as time progresses. 

Our algorithm will detect the possibility of such an event before the criticality occurs, when ac 
ceases to be a-stable (or even before this happens). It will then remove this edge from the stable 
subgraph, so the actual cocircularity will not be recorded. The new edge bd will then be detected 
by the algorithm only when it becomes at least a-stable (if this happens at all), and will then enter 
the stable Delaunay graph. In short, critical cocircularities do not arise at all in our scheme. 

As noted in the introduction, a Delaunay edge ab (interior to the hull) transitions from being 
a-stable to not being stable, or vice-versa, when the sum of the opposite angles in its two adjacent 
Delaunay triangles is n — a (see Figure [2]). This shows that changes in the stable Delaunay graph 
occur when the "cocircularity defect" of a nearly cocircular quadruple (i.e., the difference between 
7T and the sum of opposite angles in the quadrilateral spanned by the quadruple) is between a 
and 8a. Note that a degenerate case of cocircularity is a collinearity on the convex hull. Such 
collinearities also do not show up in the stable Delaunay graph. A hull collinearity between three 
nodes a, b, c is detected before it happens, when (or before) the corresponding Delaunay edge is no 
longer a-stable, in which case the angle Zacb, where c is the middle point of the (near-)collinearity 
becomes tc — a (see Figure [T3]a)). Therefore a hull edge is removed from the SDG if the Delaunay 
triangle is almost flat. The edge (or any new edge about to replace it) re-appears in the SDG when 
it becomes ca-stable, for some 1 < c < 8. 










b 


n — a 


(a) 



Figure 13. (a) The near collinearity that corresponds to a Delaunay edge ceasing to be a-stable. (b) A set of points for 
which the number of a-stable edges in DT(P) (those corresponding to the vertical Voronoi edges) is close to n. 
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SDGs are not too sparse. Let P be a set of n points in the plane. We give a lower bound on 
the number of a-stable Delaunay edges in the Delaunay triangulation of P. Our lower bound 
approaches n as a decreases to zero. 

Let «o be the number of points with no incident a-stable edges in DT(P) and let n\ be the 
number of points with a single incident a-stable edge in DT(P). Clearly the total number of a-stable 
edges in DT(P) is at least 


2(h - H 0 - Hi) + Hi _ /2H 0 + Hi \ 

2 ) 


(3) 


We now derive an upper bound on 2n °£ ni . Consider a vertex v with no incident a-stable edges. 
If v is not a vertex of the convex hull then its degree in DT(P) must be at least 27i/a (the boundary 
of its cell in VD(P) contains at least 2n /a a-short edges). If v is a vertex of the convex hull then 
its degree must be at least ( 2tc — d(v))/cc where d(v) is the angle between the two infinite rays 
bounding the Voronoi cell of v. Similarly, consider a vertex v with one incident a-stable edge. If 
v is not a vertex of the convex hull then its degree must be at least zr/a + 1 and if v is a vertex of 
the hull then its degree is at least (tt — d(v)) / a + 1. Since £d(zz) over all hull vertices is 2zr, we get 
that the sum of the degrees of the vertices in DT(P) is at least 


n i 


n 

- +1 

a 


2 H 0 


71 

a 


2tc 

a 


(4) 


On the other hand, the Delaunay triangulation of any set with h vertices on the convex hull consists 
of 3h — h — 3 edges so the sum of the degrees is 6 n — 2h — 6. Combining this with the lower bound 
in Q we get that 

/ 7i \ _ tc 2n „ „ 

Hi (—hi) +2 h 0 -< 6H, 

\ a / a a 


which implies that 


2Hq + Hi \ 6H 


zr/a 


2 . 


Substituting this upper bound in Equation |3j we get that the number of a-stable edges in DT(P) is 
at least 



-2 . 


This is nearly tight, since, for any a, there exist sets of n points for which the number of a-stable 
edges is roughly h; see Figure [l3|b). 


Closest pairs, crusts, ft- skeleta, and the SDG. Let /3 > 1, and let P be a set of n points in the 
plane. The j6 -skeleton of P is a graph on P that consists of all the edges pi] such that the union of 
the two disks of radius (/3/2)d(p, q), touching p and q, does not contain any point of P \ {p,q}- 
See, e.g., BUB] for properties of the /3-skeleton, and for its applications in surface reconstruction. 
We claim that the edges of the /3-skeleton are a-stable in DT(P), provided /3 > 1 + Q(a 2 ). Indeed, 
let pq be an edge of the /3-skeleton of P, for /3 > 1. Let Ci and C 2 be the centers of the two empty 
disks of radius (/3/2 )d{p,q) touching p and q; see Figure [l4|a). Clearly Zcipq = Zc 2 pq. Denote 
6 = Zc\pq = Zc 2 pq. Each of p,q sees the Voronoi edge e vq at an angle at least 29, so it is 20-stable. 
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(a) 


(b) 


(c) 


Figure 14. (a) The skeleton edge pq is a stable edge, (b) An edge ab of the Relative Neighborhood Graph that is not 
stable, (c) Point p is disconnected in the SDG. VD is drawn with dashed lines, SDG with solid lines, and the remaining 
DT edges with dotted lines. 


We have 1//3 = cos 0 ~ 1 — 0 2 /2 or (5 = 1 + 0(0 2 ). That is, for (5 > 1 + Q(a 2 ), with an appropriate 
(small) constant of proportionality, pq is a-stable. 

A similar argument shows that the stable Delaunay graph contains the closest pair in P(f) 
as well as the crust of a set of points sampled sufficiently densely along a 1-dimensional curve 
(see (HO for the definition of crusts and their applications in surface reconstruction). We sketch 
the argument for closest pairs: If ( p,q) is a closest pair then pq G DT(P), and the two adjacent 
Delaunay triangles A pqr + , A pqr~ are such that their angles at r + , r~ are at most n /3 each, so e pq 
is (7r/3)-long, ensuring that pq belongs to any stable subgraph for a sufficiently small; see 0 for 
more details. We omit the proof for crusts, which is fairly straightforward too. 


Remark. Stable Delaunay graphs need not contain all the edges of several other important 
subgraphs of the Delaunay triangulation, including the Euclidean minimum spanning tree, the 
Gabriel graph, the relative neighborhood graph, and the all-nearest-neighbors graph. An illustration 
for the relative neighborhood graph (RNG) is given in Figure 14 (b). Recall that an edge pq is in 
RNG if there is no point r G P such that max{||pr||, ||^r||} < \\pq\\- As shown in figure that pq 
is an edge of RNG, but the angular extent of the dual Voronoi edge e pq can be arbitrarily small. 
As a matter of fact, the stable Delaunay graph need not even be connected, as is illustrated in 
Figure [l4|c). 


6 Conclusion 

In this paper we introduced the notion of a stable Delaunay graph (SDG), a large subgraph of the 
Delaunay triangulation, which retains several useful properties of the full Delaunay triangulation. 
We proved that a 4a-stable edge in (the Euclidean) DT(P) is a-stable in DT^(P), where Q is a 
regular /c-gon for any k > 2zr/a and a is the stability parameter, and that the dual Voronoi edge 
e pq in VD <> (P) contains at least eleven breakpoints. Using these properties and the kinetic data 
structure for DT <> (P) developed in the companion paper (3|, we presented a linear-size KDS for 
maintaining a Euclidean (8a, a)-SDG as the input points move. The KDS processes only a nearly 
quadratic number of events if the points move along algebraic trajectories of bounded degree, 
and each event can be processed in 0(log n) time. We also showed that if an edge is stable in the 
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Delaunay triangulation under the Euclidean norm, it is also stable in the Delaunay triangulation 
under any convex distance function sufficiently close to the Euclidean norm, and vice versa. 

Proving a subcubic upper bound on the number of topological changes in the Euclidean 
Delaunay triangulation for a set of moving points still remains elusive (in spite of the recent 
progress in Ifl9l[20ll ). but our result implies that if the true bound is really close to cubic, or just 
significantly super-quadratic, then the overwhelming majority of these changes involve edges 
appearing and disappearing while the four vertices of the two triangles adjacent to each such edge 
remain nearly cocircular throughout the entire time in which the edge exists. 

We conclude by mentioning two open problems: 

(i) Is there a KDS for maintaining a triangulation of the entire convex hull of a set of moving 
points in the plane, which is an approximate Delaunay triangulation, defined appropriately, 
and which processes only a near-quadratic number of events? In particular, can the SDG 
maintained by our KDS be extended to a triangulation scheme of conv(P) (recall Figure [TJb)), 
e.g., using the ideas from the kinetic triangulation schemes presented in 1151 1151 . which also 
undergoes only a near-quadratic number of topological changes during the motion? Perhaps 
a deeper analysis of the structure of the "holes" in the stable sub-diagram may yield a solution 
to this problem, using the fact that for every missing edge, the two incident triangles form 
a near-cocircularity in the diagram. This may lead to a scheme that fills in the holes by 
near-Delaunay edges and that has the desired properties. 

(ii) What are the other large and interesting subgraphs of DT(P) that undergo only a near¬ 
quadratic number of topological changes under a motion of the points of P of the above 
kind? For instance, can one prove that there are only a near-quadratic number of changes in 
the a-shape or the relative neighborhood graph of P if the points of P move along algebraic 
trajectories of bounded degree. 
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